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Abstract 

Manpower  problems  of  an  organization  have  long  been  of  great  con¬ 
cern  not  only  for  prediction  of  the  future  personnel  inventory  at 
different  ranks  but  also  for  analysing  the  interaction  between  states 
of  ranks  and  indicating  chances  for  promotion.  This  study  is  intended 
to  construct  and  investigate  reasonable  and  tractable  models  for  man¬ 
power.  That  is,  models  that  describe  the  dynamics  of  recruitment, 
advancement,  and  separation  of  individuals  in  a  rank-structured 
system.  It  is  tempting  to  devise  approximate  diffusion  models  for 
such  problems,  in  order  to  obtain  simple  analytical  mathematical  ex¬ 
pressions  for  quantities  of  interest. 
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1.  Introduction 


Manpower  problems  of  an  organization  have  long  been  of 
great  concern  recently,  not  only  for  prediction  of  the 
future  personnel  inventory  at  different  ranks  but  also  for 
analysing  the  interaction  between  states  of  ranks  and 
indicating  chances  for  promotion.  Virtually  every 
organization  must  plan  their  resources  for  future  demands 
and  uncertainties  and  manpower  is  most  pertinent  among  these 
resources . 

Several  approaches  have  been  developed  for  analyzing 
this  kind  of  problem.  among  these,  Markovian,  or 
seni-Markovian  schemes  (1)  (2)  have  been  used  as  well  as 
mathematical  expressions  for  steady  state  values. 

In  the  real  world,  the  numbers  in  each  rank  are 
correlated  rather  than  independent,  and  so  does  between 
different  positions  in  a  rank,  e.g.  there  night  be  a 
relatively  favored  type  of  individual  as  well  as  unfavored 
ones. 


This  study  is  intended  to  construct  and  investigate 
resonable  and  tractable  models  for  manpower.  That  is  , 
models  that  describe  the  dynamics  of  recruitment, 
advancement,  and  separation (retirement ,  failure  to 
re-enlist,  or  discharge)  of  individuals  in  a  rank-structured 
system. 

It  is  tempting  to  device  approximate  diffusion  models 
for  such  problems,  in  order  to  obtain  simple  analytical 
mathematical  expressions  for  quantities  of  interest. 

In  this  thesis  four  stochastic  models  which  describe 
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the  manpower  stocks  or  inventories  at  each  rank  are 
developed.  The  simplest  and  most  basic  model  is  based  on 
assumed  state  independence  and  is  presented  in  section  2A. 
In  the  second  model,  described  in  section  2B,  these  two 
assumption  are  relaxed.  More  realistic  models  that  allow 
for  preferencial  weighting  factors  are  presented  in  section 
2C  and  2D. 

For  each  case,  the  long  term  prediction  model  is 
developed  evolution  of  the  process  {Q ( t) ,  t>0}  ,  following 
the  procedure  of  analytically. 

Finally  we  present  the  conclusion  in  chapter  3  along 
with  some  recommendations  for  further  work  on  this  topic. 
Appendices  are  presented  for  the  algebraic  derivations. 
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2. 


Models. 


A.  Model  I 


In  a  closed  system,  recruitment  is  allowed  only  for  the 
lowest  rank,  and  advancements  are  considered  for  the  very 
next  higher  rank,  but  separation  can  take  place  from  all 
ranks. 

In  this  model,  the  organization  has  a  sufficient  amount 
of  manpower  resources,  and  the  capacity  (ceiling)  of  each 
rank  is  unlimited. 

The  model  structure  is  depicted  in  Figure  1. 


Fig  1. 

a  (t)  is  the  recruitment  rate  for  the  first  grade. 

The  personnel  who  are  in  the  rank-1  ,Q,(t)  ,  can  be 
advanced  to  the  next  grade  by  the  rate  of  u,  (t)  or  separated 
from  the  system  by  the  rate  of  v(  (t)  and  so  on  for  the 
succeeding  grade. 

The  manpower  stocks  at  time  t  which  are  in  grade  i, 

Q^(t)  are  advanced  with  probability  u^(t),  and  separated  with 
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probability  v.  (t)  ,  over  the  interval  (t,t+dt)  . 

We  will  write  down  stochastic  differential  equations 
of  the  same  form  as  that  given  by  Gaver  and  Lehoctlcy.  (3) 


Notation 

i 

dQ.(t)  =Q.(t+dt) -Q.(t)  is  used, 
i  *  * 

For  grade  1  : 


^ :  /  V  r'- 

-t 


dQ,  <t)  =Q.a(t)  dt-v,  (t)Q,  (t)dt-u,  (t)Q,  (t)  dt 
+ /qTT t)  dW*  (t)  -  JQ,  (t)  v,  (t)  dWV|  (t) 

-JQ,  (t)  ut  (t)dHU|(t) 

.  (1,1) 


For  grade  i,  i=2,3,...,!t 


dQ.  (t)=uu(t)Q..1{t)dt-  {u.  (t)+v.  (t)  }  Q ;  (t)  dt 

+Ju;Jt)Qi-i(t)  dwaH(t)  "  /u;  W  Q  ;  (fc)  dwu.(t> 

-  Jv.  ( t )  Q .  (t)  dWy  (t) 

'  .  (1,2) 


Xn  equations  (1,1)  and  (1,2)  the  terms  dWA(t)  ,  d.f^(t)  , 
dw0(t)r  dW7;(t),  and  dW0.(t)  are  the  derivatives  of 
independent  standard  Wiener  precesses,  and  such  as  are 
usually  called  Gaussian  White  Noise. 


The  increment  of  the  population  in  the  1— st  gi.ade  is 
recruitments  minus  retirements  and  promotions  in  (t , 
t+dt) .  For  small  dt  the  fluctuating  (  diffusion  ) 

component  of  recruitment  is  represented  by  a  (  '-)  ^ ' 

where  the  scale  factor  is  the  standard  deviation  of  a 
Poisson  process  with  mean  Q#a(t)  and  dW^(t)  is  N(0,dt.)  and 
the  other  noise  terms  due  to  promotion  and  retirement  are 
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obtained  from  the  same  rationale. 


To  study  the  random  noise  component,  we  introduce  the 
standardized  noise  variables  : 

X;  <t)={Q.  (t)-Q,  g.  (t)  }  /  Sty,  i=1,2,...,k 
Qitl-  ^  (1/3) 

where^q,  (t)  is  the  mean  valuenof  Q;  (t)  j  s-h>;  />_  ■ 

The  equations  (1,1)  and  (1,2)  substituted  by  (1,3)  can 
be  arranged  as  a  system  of  k  ordinary  differential  equations 
for  the  mean  value  approximation,  and  a  system  of  n 
stochastic  differential  equations  for  the  noise  components  : 


For  the  mean  value  approximation, 

tl)  -  -i  ^t)  = 

dq  |  (t)  /dt=a  (t)  -  (v ,  (t)  +u,  (t)}c(  (t) 


dq;  (t)  /dt=u;.,(t)  q;.,(t)  -  {U  ;  (t)+v;  (t )  }  q .  (t) 

1-2 , 3 , . .  .  ,n 


(1/U) 


For  the  noise  component  , 

dX,  (t)  =-  {v,  (t)  +u  ,  (t)  }  X,  (t)  dt+/a  (t>  dW*  (t) 

**/▼,  (t)q,  (t)dWVi  (t)  -/u,  (t)  q,  (t)  d Hu  (t ) 


dX;  (t)  =U.M(t)  X[.,(t)  dt-  {U  ;  (t)  +V;  (t)}X.(t)  dt 

+/U;_.<t)qiH(t)  dW„M(t)  -/*.  (t)  q.  (t)  dwu.  (t) 
-/V;  (t)  q.  (t)  dW,.  (t) 

.  (1/5) 

We  express  (1,5)  as  a  simple  form  ; 
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dX,  (t)  =-b,  (t)  X,  (t)  dt+y/a  (t)  +b ,  (t)  g  (  (t)dw,  (t) 


dX;  (t)=u;.i(t)X..,(t)dt-bi.  ( t )  X*  (t)dt 


+b.  (t)  q.  (t)dW.  (t) 


i=2,3, . . 


(1,6) 


where  b.  (t)=u.  (t)+v.  (t)  and  dW  j  (t)  is  the  standard  Wiener 
process  associated  with  grade  i,  i=1,2,... 


(Derivation  :  Appendix  A) 
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jjg  thsn  writs  noiss  variablss  in  v  sc  tor  fashion  as  , 


x(t)  =  {xt  (t) ,  x^(t)  ,  ...,xn(t)}' 

W  ( t )  —  { M  (  (t )  ,  W  x  ( t)  ,  •  *  -  #^rv(t)  }  * 


4Xli)  = 


-b,U), 
a, it)  , 
o 


-  b4it)  , 


U»W  ,  -t.lt) 


.  0  ,  ,  b„U) 


X(t)4t 


0  + 


--  ^Wi-W’KUt) 


d  W(t) 


/ 


(1/7) 


For  the  mean  value  of  each  grade  we  can  solve  equation 
(1,4)  grade  by  the  grade  or  by  the  computer^ 

i.  e. 


q  (t)=  (exp  {-^  (v,  (x)  +  u  ,  (x)  )  dx}  a  (s)  ds 

*  t 

+  q,  (0)  exp  {“ J( v ,  (x)  +u  ,  (x)  )  }  dx 

0 

g  (t)  =  Jexp{-f(v,  (x)  +u,  ( x)  )  dx }  u,  (s)  q  ,  (s)  ds 

•  i  t 

+  q^  (0)  exp  {-J(ut  (x)  +vt  (  x)  )  dx} 


(1,8) 


For  the  noise  component  we  can  find  X,  (t)  from  (1,5) 


l.  e. 


X,  (t)  =  ^ exp  {*(( v ,  (x)  +  u  ,(x)  )  dx} 

♦  s 

•  >/a  (s)  +  (v,  (s)  +u  x  (s) }  q(  (s)  dW,  (s) 

+X,  (0)  exp  (-J(v  ,  (x)  +u  ,  (x) )  dx} 

*  .  (1,9) 

and 

2{X,  (t)  }=X  ,  (0)  exp(-f(v,  (x) +u,  (x)  )  dx} 

O 

Var  {X  ,  (t)}=^exp{-2(  (v,  (x)  +u,  (x)  )  dx} 

»  t 

•  (a  (s)  +  (v,  (s)  +u(  (s)  )  q  (  (s)  }  ds 

For  the  special  case,  where,  a(t)=a,  v.(t)=v,  and  u.(t)=u, 
after  a  lcng  period  of  time (t4«)  ,  these  values  are: 


S{X,  (t)  }  =0 

Var  (X,  (t) }  =a/  (u  +  v) 
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From  (1,8)  we  can  observe  the  mean  number  of  people  in 
1-st  grade  as  t4e©  is  a/(u,+v,).  These  results  are 

analogous  to  an  tt/M/co  queueing  model  in  which  the  mean 
number  of  people  in  the  system  is  a/(u,  +v,  )  and  the 
probability  that  there  are  n  people  is  Poisson  (a/ (u,  +v,  ))  . 

(Ref.  (4)  ) 

For  subsequent  grades  we  have  to  consider  the 
correlation  between  X-  (t)  ,  X.  (t)  ,  (i#j)  .  Since  we  know 

that  these  noise  terms  are  sums  of  Normal  random  variables 
whose  mean  is  zero,  it  follows  that  3  (X  (t) } =0  and  the 
variance  is  equal  to  the  second  moment.  We  describe  the 

moment  functions  as  : 

m,  (t)  =  E  (X*(t)  } 

mx  <t)=2{Xf(t)  } 

m|1(t)  =  E{X,  (t)  Xt  (t)  }  etc. 

and 


{m.  (t+dt)  -m.  (t)}/dt  =  dm;  (t)/dt  as  dt  =  0 

i=1 ,  2  , . . . , k 

Substitute  X  j  (t+dt)=X.  (t) +dXj  (t)  into  (1,5),  square  both 

sides,  and  take  the  expectation,  than  we  can  get  moment 
functions. 

In  this  procedure,  2[  [</a  (t)  dW^  (t)  }  2  ]  is  a(t)dt, 
Var[dW^(t)  ]=E[  (dW*  (t)  }  2  ]=dt  since 
dH«.(t)  —  N  (G,  dt)  . 
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For  grade-1  and  grade-2, 


v/  n/  v/ 

dm,  (t)/dt=a(t)  +  v,  (t)  q,  (t)+u,  (t)  g  ,  (t) 

-2{v,  (t)  +u ,  (t)}a  ,  (t) 


dmz  (t)  /dt-=u  ,  (t)  g  ,  (t)  +  (t)  +  v4  (t)  }gx  (t) 

-2  {ut  (t)  (t)  (t)  +2ut  (t)  m(Jt) 


dm|t(t)  /dt=u  ,  (t)  m,  (t)  -  {u  l  (t)  (t)  +u  ,  (t)  +v,  (t)  }  m(1(t) 


-u  ,  (t)  g,  (t) 


(1,10) 


(Derivation  :  Appendix  3  ) 

The  variance  and  covariance  are  the  solutions  of  these 
equations.  These  solutions  can  be  determined  numerically, 
on  computer,  and  steady  state  values  can  be  derived 
explicitly . 
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The  graphical  expressions  are  available.  (  graph-1) 


THE  STEADY  STATE  MEAN  IN  EACH  RANK 

OF 

A  MULTI-RANK  ORGANIZATION 

When  an  organization  maintains  a  certain  policy  for 
recruitment,  promotion,  and  retirement  indefinitely,  it  may 
have  a  steady  state  property  (when  the  rate  of  flowing  into 
a  grade  is  less  than  that  of  flowing  out) ,  thus  an 
analytical  approximation  technique  can  provide  an  answer 
that  is  quite  simple  and  useful. 

For  the  purpose  of  investigating  such  a  situation  we 
take  a  (t)  ,  Uj  (t)  ,  v.  (t)  to  have  constant  values  in  equations 
(1,4),  and  solve  these  under  steady-state 

conditions.  This  can  be  done  by  setting  the  derivatives 

equal  to  zero,  and  solving. 

q ,  (»)  =a/b, 

q^  («o)=u,  a/(b,  \  ) 

•  •  • 

<3h  (<o)  =u,  Uj . . .  u^a/fb,  bt  .  . .  bn } 

where  b;  =u.  +v.  ,  i=1,2,...,n 
Special  Case 

If  there  is  the  same  chance  of  promotion  and  retirement 
for  every  grade,  u  and  v,  then 

g  (»)  =  (a/u)  {u/b}1,  n=1,2,... 
n 

and  the  total  number  of  people  in  the  system  is 

2  g.  (»)  =  (a/b)  {[  1- (u/b)"  ]/[  1-  (u/b)  ]}  ,  where  (u/b)  <  1 

;-i  < 
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ri 

Let  g  =  X  g.  (oo  )  ,then  g  is  the  upper  level  of  total  number 

n  ;*»  1 

of  people  up  to  rank-n,and  the  maximum  number  of  people  in 
this  organization  is 

q  =  1  g.  <«)  =  (a/b)  {1-  (u/b)  }-i  ,  where  (u/b)  <  1 
«  ' 

Sometimes,  one  is  interested  in  the  number  of  people  who  are 
above  a  certain  rank. 

g  -q  =  (a/b)  ( (u/b)*  /[  1-  (u/b)  ]}  ,  a=1,2, .  .  . 

fO  n 

.  (Mi) 

This  geometrically  decreasing  function  represents  the 
relative  difficulty  in  occupying  the  higher  ranks.  Moreover, 
the  geometric  growth  rate  (u/b)  might  be  of  interest  in 
comparing  organizations. 
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B.  Model  II 


In  this  model  we  have  improved  a  ceiling  for  the 
population  size  of  each  grade,  it  is  not  allowed  to  exceed 
that  level,  and  promotion  rate  is  proportional  to  the 
vacancy  of  the  very  next  higher  grade,  where  the  vacancy  is 
the  difference  between  ceiling  and  present  population. 

Let  u  j  (t)  be  the  rate  at  which  individuals  become 
candidates  for  promotion  from  grade  i.  Some  of  those 
in  a  grade  can  be  promoted,  and  the  rest  of  them  will  "give 
it  up",  i.e.  be  forced  to  leave.  We  have  also  a  natural 
loss  from  each  grade,  described  by  the  rate  v.  (t)  . 

The  model  structure  is  depicted  in  Fig  2. 

Q  :  ceiling  of  i-th  grade 


Y 


Fig-2 
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The  probability  of  being  recruited  from  the  population 
is  a  (t)  ( t ) /QJ} ,  the  probability  of  being  promoted  is 
modeled  by  u  (t)  (1-[QX  (t)  /Qj}  ,  and  the  probability  that 
present  grade  occupants  give  up  this  organization  because 
they  failed  to  advance  as  candidate  is  u,  (t)  {(^(lO/Q  }  . 

In  order  to  express  the  change  in  a  grade  we  shall 
first  study  the  change  in  the  mean.  The  mean  increment 
to  the  grade  is  the  recruitment  or  advancement  into  it., 
minus  those  separated  and  becoming  eligible  for  promotion  in 
(t,  t  +  dt)  . 

Concerning  fluctuation  (  diffusion  ) ,  we  have  to 
consider  the  terms  due  to  recruitment,  promotion,  giving  up 
the  organization,  and  natural  loss  separately  since  some  of 
these  terms  are  correlated  with  ones  in  the  next  grade. 

The  noise  term  dWg^ 1 ^denotes  the  fluctuation  due  to 
losses  from  the  organization,  and  dWu  t 2  \t)denotes  that  due  to 
promotion  to  the  next  grade. 


17 


For  the  1-st  grade 


dQ,  (t)  =  Q0  a  (t)  {1-[Q,  (t)/Q,]}dt-u,  (t)Q,  (t)  {Qi  (t)/Q;dt 
-u,  (t)Q,  (t)  { 1 — C  Q a  (tj/QJJdt-v,  (t)  Q,  (t)  dt 
+/Q.a(t)  {1-CQ,  ( t )  / Q,J}  dWA  (t) 

Vu,  (t)  Q,  (t)  {Q,  (t)/Q,  }  d » <  i  >  (t ) 

-/u,  (t)  Q,  (t)  (1-CQ,  (t)/Qj}  d  2  >  (t) 

-/v,  (t)  Q,  (t)  dWV|  (t) 


where  Q.  (t)  <Q.  ,  i=1, 2, . . . 


<2, 1) 


For  the  2-nd  grade. 


dQ,  (t)  =u  ,  (t)  Q  ,  (t)  {1-[Q,  <t)/Q,]}dt 


-ua  (t)  Q,  (t)  {Q 3  (t)/Q,  }dt 


-u,  (t)  Q,  (t)  (1-C  Q3  (t)  /Qj]}  dt-v,  (t)  Q,  (t)  dt 


Vu,  (t)Q,  (t)  {1-[Q,  (t)/QJ)  d tf  <  2  >  (t) 


(t)  Q,  (t)  {Q  3  (t)/Q3J  d»<  *  >  (t) 

2 

■/ul(t)Qx  (t)  c  1-C  Q3  (t)/Q,3)  d^C  2)(t) 


-A  (t,  Q,  (t)  d tf Vi ( t ) 


where  Q,  (t)  <Q,  ,  Q  3  (t)  <Q, 


(2,2) 
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We  standardize  noise  variables  and  Q  as  k  Q  to  apply 
the  same  technique  which  we  had  before. 

(Appendix  A  ) 

We  can  get  two  deterministic  mean  value  functions  and 
two  stochastic  differential  equations  ; 


dq,  (t)/dt=a(t)-{a(t)/k,  +u,  (t) +v,  (t)  }  q  ,  (t) 


dq,  (t)/dt=u  ,  (t)  g,  (t)  -  {  u  t  (t)  +  v*  (t)  }qt  (t) 
,  (t)q,  (t)g  ,  (t)  A» 


.  (2/3) 

and 


dX  ,  (t)  =-  {u  ,  (t)  +  v ,  (t)  +[  a  (t)  A  ,  ]}  X  ,  (t)  dt 
+/a(t)  {1-[q,  (t)/k,  ]}  dW*  (t) 

-/g,  (t)  q;  (t) u ,  (t)/kt  dWO ) 

-/u,  (t)q,  (t)  (1-Eq^t)  A*  ])  dW<  2  >  (t) 

-A,  (t)q,  (t)  dWV(  (t) 

dX2  <t)  =  {1-[gi  (t)/kt  ]}u,  (t)X,  (t)  dt 
-Cu»(t)+va  (t)}X4  (t)dt 

(t)g,  (t)  {1-[q,  (t)  Aa  B  2  >  <*) 

-/ul(t)qi(t)q3(t)/k3  dW^X  )  (t) 

*/“x  (t)ql  (t)  (i-[g,(t)  a3  3)  (t) 

~/vx(t)qt  (t)  dWVi(t) 

.  (2/4) 
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We  express  tbese  in  a  vector  fashion 

2  (t)  =  (x,  (t)  ,  x2  (t)  }  • 

E  (t)  =  £Wa(t)  ,  W<D(t),  Wtt<2>(t),  HU)  (t)  , 
W«<t2>  (t)  »  W^(t)  }• 

and 


4  X(.+)  = 


S 


{ VM  +  ui(t)| 


Xtt)  dt 


o  0  /u,wy 


0  .  ,  0 


N 


dwltj 


(2,5) 
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Since  we  are  interested  in  the  variance  of  the  noise 
factor  for  each  grade,  we  make  second  moment  functions  from 
(2,4)  including  covariance  between  X,  (t)  and  X2  (t) 

(  Appendix  3  ) 

dm,  (t)/dt=  a(t)  +  {u,  (t)  +v,  (t)}q,  (t) 

-£2a(t)q,  (t)/k,  }-2£u,  (t)+v,  (t)}m,  (t) 

dm j  (t)/dt={ua  (t)  +va  (t)  +u,  ( t )  [ q ,  (t)/kt  ]}  {gx(t)  -2m.,  (t) } 

+  U>  (t)q,  (t)  {  1**&q  x  (t)  /k  t  ]} 

+  2u,  (t)  £1-[q4(t)/k4  ]}m(x(t) 

dm11(t)/dt=u,  (t)  {l-tgjtj/k,]}  [m,  ( t )  —  g  ,  (t) } 

-£u,  (t)+ul(t)  +v,  (t)  +v2  (t) 

+  u,  (t)[g,  (t)  /k2]}  m.Jt) 

.  (2/6) 

We  can  solve  (2,6)  together  with  (2,3)  as  five 
first-order  differential  equations. 
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The  graphical  expressions  are  avaiable.  (graph-2) 


TH2  STEADY-STATE  MEAN  IN  EACH  HANK 

OF 

A  MDLTI  HANK  ORGANIZATION 


For  long  term  prediction,  an  analytic  approximate 
technique  is  very  attractive.  In  this  case  the  organization 
maintains  the  same  policy  of  recruitment  and  promotions  . 

From  (2,3) 


K«)  -  Q.  {— - } 

Q,  (  <J,  +  v, )  +  d  J 

a  Q;[ 


Q»-, 


Q;  C  U|  +  Q 
Q;-,  U._, 


^,’-iCao) 

5,.. 


3 


1*2.  3. 


We  introduce 


O: 


*«-■ 


+  V. 


L  „  6.  (  Uj  +  V.  ) 

wWe  ^  -■  1  _ '  1 

Qi-»  •d’.-i 


i  *  2  ,  3 
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and  let 


ss-  l  n 


-i 


■'  +  5;s;- , 


i  3  2  .  3, 


(2/7) 


This  formula  describes  the  mean  number  in  rank  i  as  a 
function  of  policy  and  behavioral  parameters. 


For  the  same  ceiling  of  each  grade  and  the  same  chances 
of  promotion  and  retirement,  it  is  very  simple  to  express 
the  general  term  of  the  number  of  people  in  a  grade. 

por  ^  =  $ 

j.„  U-S)  g'~'  £,*(■-  s’") 

*'  I  -  s 

i 

-  (  »  -  $  ) 

1  (~)  *  Q - : - - - - - 

S',  - 

-S.nce  §  >  i  j  -Jy-or  +He  higher  grade, 

2. 

? 


V«0  =  Q 


=  K§ 


1  +  y,  -  ? 


r 


-i 


i  *  ^  .  3 


✓  y 


where  ^  _ 


Q0  +  zr) 

i  + 


(2,3) 


2U 


This  looks  like  a  geometric 
3.  <=  3  the  distribution  of 

geometric  in  the  steady  state. 


progression.  Thus  when 
individuals  over  ranks  is 


In  the  following  section  we  develop  this  model  to 
analyse  the  interaction  between  the  two  kinds  of  individuals 
in  a  grade. 
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C.  Model  III  (Initial  preference  Model  ) 

We  modify  Model  II  for  two  kinds  of  populations,  called 
favored,  and  unfavored  (  or  normal  ) .  In  this 
organization  the  recruitment  rate  of  the  favored  population 
is  proportional  to  the  vacancy  i.e.  number  of  openings  in 
the  grade  to  which  they  are  candidates  for  advancement,  and 
for  the  unfavored  group  it  is  proportional  to  the  vacancy, 
but  also  to  the  number  of  the  members  of  the  favored  group 
who  are  employed  therein.  It  is  very  interesting  to 
observe  one  of  these  examples  in  an  organization  which  has 
male  and  female  components  ;  perhaps  the  female  is  favored 
as  a  policy  choice. 

Let  G ( t)  ,  and  the  subscript  g  ,  denotes  the  favored 
individuals,  and  H(t)  and  the  subscript  h  denotes  the 
unfavored  ones. 

The  model  structure  is  depicted  in  figure  3. 


Fig-3 
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For  favored  individuals  in  the  1-st  grade. 


dG,  (t)=a^(t)Q0  {1-[G,  (t)+H  ,  (t)  ]/Q,  }dt 
(t)G,  (t)  dt— U (t)  G  ,  (t)dt 
+l/av(t)  {1-f  G  ,  (t)  +H ,  (t)  ]/a,  }  dW^t) 


(t )  G ,  (t)  awVi<t)  -\Ai% Ttj  G,  (t)  dWUi(t) 


where  G,  (t)  +H,  (t)  <  Q, 


(3,  1) 


For  unfavored  individuals  in  the  1-st  grade, 

dH,  (t)  =ah  (t)Q0{G,  <t)/[G,  (t )  +  H  ,  (t)  ]} 

(t)+H,  (t)  ]/Q,  }  dt 
-  {vh  (t)  +uh  (t) }  H  .  (t)dt 
+/aK(t)Q9  {G,  (t)/[G,  (t)+H,  (t)  ]}" 

•  {1-[G,  (t)  +H  ,  (t)  ]/QJ  dW  (t) 

-^{tJH,  (t)dHVh(t) 

-/uh  (t)  H,  (t)dW„  <t) 

1  hl 

.  (3,2) 
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To  study  the  random  noise  component,  we  introduce  the 
standardized  noise  variables. 


For  favored  individuals, 

X.  <t)  =  {G.  (t)-Q  q.  (t)  }/75  ,  i=  1  /  2  , .  .  .  ,n 

I  I  0  |  • 

For  unfavored  individuals, 

Y.  (t)  =  £H  .  (t)  -Qo  h;  (t  )}//Qo,  i=1,2,...,n 
2;  =ki  Qo 

.  (3,3) 

where  a.  (t)  and  h.  (t)  are  the  mean  value  of  G.  (t)  and  H.  (r)  . 


From  equations  (3,1)  ana  (3,2)  substituted  by  (3,3),  we 
can  get  a  system  cf  two  ordinary  differential  equations  for 
the  mean  value  approximations  and  another  system  of  two 
stochastic  differential  equations  for  noise  components; 
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ag,  (t)/dt=aj  (t)  £1-C  g,  (t)  +h  ,  (t)  ]/k,  } 
~{Vj  (t)  +u  j.  (t)  }  g(  (t) 


dh, 


where  g( 


(t)/dt=[1/[g(  (t)+h,  (t)  ]-  1/Jc  , }  (t)g,  (t) 

~{vh  (t)  +uK  (t)}h,  (t) 

(t)  +  h,  (t)  <  k, 

.  (3,4) 


dX,  (t)  =-{a?(t)/k,}  {X,  (t)  +Yt  (t)}dt 
-  {Vj.  (t)  +u-  (t)  }  X  ,  (t)dt 
Va^(t)  [1-[g,  <t)  +  h,  (t)  ]/k, T  dWA^(t) 
Vv.  (t)g,  (t)dw^  (t)  -/u}  (t)  g,  (t)  dwu  (t) 

M  T i  fj 


dYv  (t)  =  {1/C  9 1  (t)  +h  ,  (t)  ]-1/lt,  }ah  (t)  X,  (t)  dt 

(t)dt 

-[ah  (t)g,  (t)/(g,  (t)+h,  (t)}2]{X,  (t)+Y,  (t)}dt 
+/au  (t)g,  (t)  { 1/[  g  ,  (t)  +h  ,  (t)  ]-1/k,}  d wa(|(  t ) 


-/vh  (t)  h,  (t)  dWv  (t)  -Ahi(t)  h,  (t)  dW„H(t) 
where  g  (t)  +h  (  (t)  <  k , 

.  (3,5) 

(  Derivation  :  Appendix  C  ) 
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We  then  express  the  noise  variables  in  vector  fashion 


Z (t) =  (X  ,  (t)  ,  Y,  (t)  }  ' 

W(t)  =  {Wa  (t)/  WA  (t)  ,  Wu  (t),  WU  (t)  ,  W„/t),  (t)  }* 

1  h  r,  3,  K 


d2(t)  = 


l*) 


+  ti) 


1 


-Ljv*) - r 

U'*1  - , 


■  \ 

-  - 


\ 


I  ?,«)  + 


Z  c-t)  <Jt 


4- 


a.tuj  i-  M 1  iL*> 


0 


\ 


0 


Alt) 


2,^V+)  A 


;)},w)  o 


o  ,  7°4,W  8,h=>  ,  0 

7u»rtH,rt)  ,  o  ,  /vTwiflii*) 


V, 


dwct) 


/-1 
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Since  X,  (t)  and  Y,  (t)  are  approximately  normal  random 
variables  with  zero  mean,  we  can  find  their  variances  by 
solving  moment  function  differential  equations. 

Let 


my  (t)=E[Xf  (t)  ] 
n^(t>=E[X,  (t)  Y,  (t)  ] 
mY((t)=E[Y2<t)  ] 

and 

(X2  (t+dt) -Xf  (t)  }/dt=dmx  (t)/dt  as  dt-*0 
and  similarily  for  m  (t)  and  mY(t). 

Ti 
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By  squaring  both  sides  of  equation  (3,5)  and  taking 
expectation,  we  get  a  system  of  three  ordinary  differential 
equations  of  moment  function  ; 

dmXi(t)/dt=a^(t)  { 1 — C  9  ,  (t)+h,  (t)  ]/k  } 

+  {v3.  (t)  +Ug.  (t)  }  g,  (t) 

-2  {v^  (t)  +  u^  (t)  +[  a$  (t)  /k ,  ]}  m  X|  (t) 

-2{a^(t)  /k,  Jm^t) 

diaY(t)/dt=as  (t)  g(  (t)  (1/[g(  (t)+h,  (t)  ]-1/k,} 

+  {v„(  (t)  +ub  (t)  }  h,  (t) 

+  2ah  (t)  {1/[  g  t  (t)+hk  (t)  ]-1A,}nv,Y(t) 

-2  (vh  (t)  +utt)(t)  }  m Y (t) 

-2ah  (t)g,  (t)/{g(  (t)  +h  ,  (t)  }  2  (ra^t)  +mY(  (t)  } 

dmXiX(t)/dt=ah  (t)  £  1/C  g  ,  ( t)  +h,  (t)  j-l/k.Jm^t) 

-ah  (t)  C  g,  (t)  /  £g ,  (t)+h,  <t)}2  JmXi(t) 

-  (v^(t) +vh  (t)  +u^(t)+uh  (t)+[a^(t)  A  ] 

+  aVl  (t)  g ,  (t)  /[  g,  (t)+h,  (t)  ]*}  rn^t) 
-(aj(t)/k,}^(t) 

.  (3,6) 
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In  order  to  get  the  variance  and  covariance  of  the 
noise  component,  we  can  solve  a  system  of  three  ordinary 
differential  equations  instead  of  two  stochastic  equations, 
(3,5)  ,  together  with  (3,4)  . 

For  higher  grades,  up  to  n-th  grade,  it  is  required  to 
set  up  a  system  of  2n  ordinary  differential  equations  for 
mean  value  approximation  and  another  system  of  2n+(^) 
ordinary  differential  equations  derived  from  a  system  of  2n 
stochastic  differential  equations  for  noise  component. 

Taking  a  six  grade  military  rank-structured 
organization  as  an  example,  we  are  required  to  solve  12 
equations  for  mean  value  approximation  and  78  equations  for 
noise  component,  which  makes  a  system  of  90  ordinary 
simultaneous  differential  equations. 

We  will  present  a  simple  numerical  example  for  the 
first  grade. 
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Graphical  expressions  are  available  .  (graph-3) 


THE  STEADY  STATE  I1EAN  IN  EACH  RANK 

OF 

A  MULTI  RANK  ORGANIZATION* 

For  the  special  case,  when  we  have  constant  parameters 
for  the  recruitment,  retirement,  and  promotion,  it  is 
attractive  to  develop  an  approximation  technique  for  the 
number  of  people  in  each  grade. 

To  observe  the  steady  state  value,  we  set  the 
derivertives  equal  to  zero  in  equation  (3,4),  and  solve  for 
the  ratio  of  the  number  of  favored  population  and  unfavored 
ones,  grade  by  grade. 

For  simplicity  we  adopt  the  following  notation  ; 


w.  =  g.  (coJ/hj  (oo) 

bi  =  {vv/cvv 

a=a3/aK,  u; 

By  solving  (3,4)  substituted  by  (3,7)  , 
w,  =  {a+/a2+iTo,  a  }  /  (2b,  ) 


(3,7) 


(3,8) 


For  the  general  grades. 


(3,9) 
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V'h*1'  (V»i )  }=bi  wi 


(3,10) 


and 


*;  ®f;  w.^Vcf;  w;.,+ 1)2-1  where  f.  =U;_,/(2b.  ) 

2,3,..  . 


(3,11) 


It  is  straight  foward  to  get  mean  value  approximation 
in  steady  state,  when  it  exist,  to  solve  (3,4)  and  (3,9) 
substituted  by  (3,8)  and  (3,11) 


9,  (  ®)  =a^/{u3  +v^  +  (a^  /k,  )  (1+1/w,  )} 


gi  (C0)=UJ.  gi->  *  («j.  gMA‘,  )  (1  +  1/w-,  ))/  i=2,3, 

.  (3,12) 


and 


h j  (°°  )  —  9 j  (®o)  /w;  ,  i— 1,2,... 

.  (3,13) 

For  the  special  case,  when  we  have  the  same  policy  of 
promotion  and  retirement  for  all  rank  ; 

v.  =  f  w.^+j/  (fw._  +  1)2-1  from  (3,10) 

It  is  very  interesting  to  observe  that  the  ratio  of 
g  (oo  )  and  h  (<o  )  has  a  tendency  to  converge  a  certain 
constant  for  the  very  high  rank.  We  can  track  the  value  by 
assuming  w;  =  w;_,  ,  then, 

uw  ( 1-1/w)  =bv 

w=u/  (u-b) 
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(3,14) 


where  u=u^/u^  b=  (u^ +v^  | /^  +  vK  )  and  u*b 


But  ,  when  we  have  the  same  policy  for  the  two  groups, 
f = 1/2  ,and 

w  =.5w  +/(.5w  +1)2-1 
which  diverges. 
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D.  Model  IV  (Modified  Initial  Preference  Model) 

This  model  is  based  on  the  same  general  concept  that 
was  used  to  develop  Model-Ill.  The  only  difference  between 
Models  III  and  17  is  in  the  preference  factor. 

We  use  G;  (t)  /Q.  instead  of  G;  (t)/{G.  (t) +9.  (t)  }  as  a 
preference  factor  ,  which  requires  only  simple  algebra  to 
derive  the  final  result  ,  at  least  in  the  steady  -state 
situation.  Accordingly,  (3,2)  is  changed  to  the  following  : 

dH,  (t)=ah  (t)Q0  (G,  (t)/Q,  }  C 1  C  G  x  ( t)  +H  ,  (t)  ]/2,}dt 
-{vh((t)+uh  (t)}H,  (t)dt 

+/ah(t )Q#{G,  (t)/Q,J  £1-[G,  (t)+H,  (t)  ]/Q,  }  d ( t ) 

S\  (t)  H,  (t)  dW^(t)  -/uK(t)H,  (t)dHdh(t) 
where  G,  (t) +H,  (t)<Q, 

.  (4,1) 

We  derive  a  system  of  two  ordinary  differential 
equations  for  the  mean  value  approximation  and  another 
system  of  two  stochastic  differential  equations  for  the 
noise  component  form  (4,1)  together  with  (3,1). 

dg,  (t)/dt=a^(t)  { 1  — £  g ,  (t)  +h,  (t)  ]/k, } 

-  (Vj  (t)  +u^(t) }  g,  (t) 

dh,  (t)/dt=-  {vh(  (t)  +uh  (t)  }  h,  (t)+av,(t)g|  (t)/k, 

-(ah  (t)/k2)  (g (  (t)+h,  (t)}g,  (t) 
where  g(  (t)  +  h(  (t)  <  k 
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dx,  (t)  =-  {a^(t)/k,}  (X  ,  (t)  +  Y ,  (t)  }dt 

-{V^(t)+Uj(t)}X4  (t)dt 

Va^(t)  {1-[g,  (t)+h,  (t)  ]/k,  }  dW^(t) 

Vvj.(t)g(  (t)dw^(t) 

”/u3-  (t)  g  ,  (t)dwu  (t) 

dY,  <t)  =  {[ah  (t)/k,  ]-2g,  (t)-h,  (t)}X,  (t)  dt 
-{uh  (t)+vh  (t)  +g,  (t) }  Y,  (t)dt 
+/{ah  (t)/Jc,  }  {1-[g,  (t)+h,  (t)  ]/k,  }g,  (t)dWAfc(t) 
h,  (t)  dH„h  (t)  Vu^dih,  (t)  aWy  (t) 


(4,3) 
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We  express  it  as  a  vector  fashion  ; 


Let 


Z(t)  =  {X,  (t)  ,  Y,  (t)  }  • 


then 


Z(t)  {W.  ,  w*  ,  w  ,  h„  ,  w0  ,  wv  ,  }  • 

<J  "  it  N  h> 


Jzw  = 


j  (  I 

~\  ^  4  U3,rt)  4  V*’ |  ; 


Ojt*) 

•ft. 


/ 


JW^Wl  ,  0  ,  ■ 


J  V»V*)  ,  0 

0  ,  JOfcWfit) 


4W(-t) 
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three 


By  the  same  logic  a  system 
differential  equations  for  moments  is 


of 

derived 


from 


dmx>(t)/dt=a^(t)  {1-[g,  (t)+h,  (t)  ]/k,  } 

+  {v^  (t)  +u^(t)  }  g(  (t) 

-2{Vj  (t)  +uv  (t)  +[a^(t)/k,  ]}  aX|(t) 

-2{a^  <t)/k,  }mM(t) 

dmxx(t)/dt=[[ah(t)/kl  ]-2g,  (t)-h,  (t)}mx  (t) 

-  {a^(t)/k,  }  {mY  (t )  +n»x^(t)  } 

-{v^(t)  +vh  (t)  +Uj  (t)+uh  (t)+g|  ( t)  }  mx (^(t ) 

dmx(t)/dt={1-[g(  (t)+h,  (t)  ]/k,  }  {ah  (t)/k,  }g,  (t) 

£vh|  (t)  +uh  (t)  }  h.  (t) 

+  2{[ah  (t)/k,  ]-2g,  <t)-h,  (t)  }  m^t) 

-2{g,  (t)  +vh  (t)  +uHi(t)}mYi(t) 

.  (4,4) 

He  solve  (4,4)  together  with  (4,2)  to  get 
value  approximation  and  variance  and  covariance  of 
component. 


ordinary 
(4,3)  . 


the  mean 
the  noise 
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NUMERICAL  EXAMPLE 
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Assuming  zero  initial  values,  then 
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4  3 


3. 


CONCLUSION 


In  this  thesis  ,  we  have  developed  four  Models  to 
analyze  man-power  stocks  and  flows  in  a  rank  -  structured 
hierarchy . 

We  applied  the  diffusion  approximation  technique  to 
investigate  the  fluctuation  due  to  recruitment,  retirement, 
and  promotion  rather  than  Markovian  or  semi-Markovian 
schemes . 


From  the  very  simplistic  Model-I,  which  has  an 
unlimited  ceiling  and  independent  transition  policy,  we 
developed  Model-II  whose  recruitment  and  promotion  policy 
are  dependent  upon  the  present  population  of  the  very  next 
grade. 


The  mean  value  approximation  and  variance  of  the 
component  (  fluctuation  )  is  the  essential  facto 
planning  purpose  as  well  as  prediction.  It  is  very 
to  analyze  the  relationship  between  the  fluctuations  o 
tandem  grades  as  the  latter  depends  upon  the  pro 
policy  (  rate  ) . 


noise 
r  for 
useful 
f  two 
motion 


Model-Ill  and  Model-IV  were  developed  to  represent  a 
more  complicated  organization  than  that  of  Model-II.  In 
this  model  two  individual  types  (  favored  and  unfavored  ) 
compete  for  the  same  positions.  The  relationship  between 
neighboring  grades  as  well  as  unfavored  and  favored 
individuals  in  the  same  grade  were  analyzed. 


The  approximation  technique  for  the  steady  state  case 
can  be  applied  to  the  real  world  easily  because  of  its 
simplicity. 

Finally,  the  man-power  problem  has  been  concerned  with 
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partial  change  (modification)  in  an  organization,  rather 
than  with  the  whole  structure.  In  order  to  analyze  these 
effects,  with  their  fluctuations  resulting  from  a 
modification  of  a  part  of  the  policies,  it  is  very 
attractive  to  develop  models  using  diffusion  approximation 
techniques.  They  are  relatively  simple  and  show  promise 
of  representing  realistic  systems  with  resonable  accuracy. 
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Appendix-A  (Derivation  ) 


Substitute  (1,3)  into  (1,1)  and  rearrange  it  for  dx,  (t)  ; 

i 

dx,  (t)=-{v*  (t)+ut  ( t) } x ,  (t)dy 

VQ.  {a(t)dt-[v,  (t)  +u,  (t)  ]q  ,  (t)  dt-dg,  (t)  } 

-/V,  (t)  {[X,(t)//Q,  ]+g  ”  (t) }  dWVi  (t) 

Vu,  (t)  {[X,  (t )//Qa  3+g,  (t)}  d wu  ( t ) 

+l/T(t7dWa(t) 

.  (A-l) 


Since  X,  (t)  is  normal  random  variable,  dX,  (t)  is  finite 
,  and  v/Qo  term  in  (A-1)  must  be  zero  as  Q# oo  (  for  a  huge 
population  ) ,  otherwise  it  will  blow  up. 

It  follows  that  ; 


dq,  (t)  /dt=a  (t)  -  {v,  (t)+u,  (t ) }  g  t  (t) 

and 


(A-2) 


dX  (t)=-{v,  (t)+u,  (t) }  X ,  (t)dt 

+/a  (t)  dW*  (t)  -A,  (t)  q  ,  (t)  dWv  (t) 

Vtt,  (t)q,  Ct7dWu  (t) 

.  (A-3) 
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For  the  general  grade  we  rearrange  (1,2)  substituted  by 
0,3)  . 


dX.  <t)=U._t(t)X;_(t)dt-{U;  (t)  +V.  (t)}X;  (t)dt 

VC.  {«^(t)qH(t)dt-CU,  (t)  +v.  (t)  ]g  .  (t)-dq.  (t)  } 
+A_,(t)  {[X.JtJ/yQ.  ]+q._t(t)  }  dWUfJt) 

-/Uj  (t)  {[X.  (t)//Q,  ]+q.  (t)  }  dHu.  (t) 

~/T(t)  ([X.  (t )//Q  ]+q.  (t)  d Wv.  ( t ) 

1  •  V  |  | 

.  (A-4) 


By  the  same  rationale  which  we  had  in  (A-1)  as  Q  — »co  ; 


dq,  (t)/dt=u._((t)q.  ((t)  -{u.  (t)  +  v.  (t)}q.  (t) 


(A-5) 


and 


dX;  <t)=U;.,(t)XU|(t)dt-CU;  (t)+V;  (t))X;(t)dt 
+A.,(t)qi.l(t)d^J;  (t) 

-/u5  (t)q.  (t)dWu.(t) 

-/V;  (t)  q.  (t)dwv.(t) 
i=2, 3, . . . 

.  (A-6) 
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Since  noise  terras  in  a  grade  are  independent  of  each 
other,  and  the  standard  Wiener  process  dW (t)  is  N(0,dt)  and 
its  scale  factor  is  the  standard  deviation  of  a  Poisson 
process,  we  can  express  (A-3)  and  (A-6)  in  a  simple  form, 

i.  e. 


v/dfdRCk(t)+v/d2dWk  (t)=/d2  +  d2dW(t) 
therefore 

dX,  (t)=-b,  (t)  X,  (t)  dtVa  tt)  +To,  (t)  g  ~  (t)dW,  (t) 

dX;  (t)=u..,(t)  X;_(t)dt-b.  (t)  X .  (t)dt 

+/u..((t)gi_i(t)+b.  (t)  g ,  (t)dW.  (t) 
i=2,3 ,  . . . 

where 


b;  (t)=v;  (t)  +  U.  (t) 
i=1,2, . . . 


(A-7) 


The  subscript  for  the  Wiener  process  represents  its 
grade. 
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Appendix  B 

(Derivation  of  Moment  Function  ) 

To  solve  the  variance  and  covariance  of  the  noise 
component,  re  derive  a  system  of  simultaneous  ordinary 
differential  equation  for  moments  from  the  stochastic 
differential  equations  (1,6). 

X;  ( t+dt)  =X .  (t)+dX;  (t) 

.  (B-1) 

Substitite  (B-1)  into  (1,6)  and  square  both  sid<2s  and  take 
the  expectation. 

1,  (t)=E[X,2(t)  ] 

m  !  (t)  =E[  X2  (t)  ] 

mli(t)=E[X,  (t)  X,  (t)  ] 

and 


{E[X,2(t+dt)  ]-E[X,2(t)  ]}/dt=dm,  (t)/dt 
as  dt  approaches  to  zero. 
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For  the  1-st  grade  ; 


and 


X,  ( t+dt )  =X ,  (t)-b,  (t)  X,  (t)  dt  +  /ut)+b,  (t)q,(t)  dW,  (t) 

.  (B— 2 ) 


X2(t+dt)=X2(t)+b2(t)X2  (t)  (dt )  2 

+  {a  (t)  +b,  (t)  g,  (t)  }dt-2b,  (t)  X,2  (t)  dt 
♦  2/a(t)+b,  (t)q,(t)X,  (t)dW,  (t) 

-2b,  (t)y/a(t)  +b  j  (t)  g  ,(t)  X,  (t)  dW,  (t) 

and 


(2[X2(t+dt)  ]-E£X2(t)  ]}/dt  : 

b,2  (t)E[X,2(t)  ]dt 

+  {a(t)+b,  (t)  q  t  (t)}-2bl  (t)E[X2(t)  ] 

+  2/a  (t)  +b  ,  (t)  q  ”  (t)  E[X,  (t)dW,  (t)  ]/dt 
-2b,  (t)/a  (t)  *b,  (t)  q  ,  (t)  E[X,  (t)dW,  (t)  ]/dt 

Since  X,  (t)  and  dW,  (t)  are  independent,  its  covariance 
is  zero,  and  "dt'*  term  is  zero  as  "dt"— >0 

dm  ,  (t)  /dt=a  (t)  +b  ,  (t)  q  ,  (t)  -2b  ,  (t)  m  ,  (t) 

.  (B— 3 ) 

-  l.C*)-  2.  (v,+ 
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By  the  same  rationale  we  derive  covariance  function. 


Vx(t)' 

From  (1,6)  substituted  by  (B-1) 

X,  (t+dt)  =X ,  (t)-b,(t)X,  (t)  dt+/a  (t)  +b  ,  (t)q(  (t)'dw  ,  (t) 

.  (B-4) 

Xt  (t+dt)  =Xx  (t)  +  u  ,  (t)  X  ,  (t)  dt-bx  (t)  X^  (t)  dt 

+/u,  (t)qt(t)+bt(t)qjt)'dw*(t) 

.  (B-5) 

Take  the  product  (B-4)  and  (B-5)  and  get  dm.jtj/dt  by 
the  same  rationale  which  we  had  in  (B—  3) 

Since 

X,  (t)  and  da ,  (t) 

X  ,  (t)  and  dW,  (t) 

X x (t)  and  da,  (t) 

X  L(t)  and  dWt  (t) 

dW,  (t)  and  da*  (t) 

are  independent 
and  E[  X  ,  (t)  dH  ,  (t)  ]=0,  etc. 

dm,  x(t)  /dt=u  ,  (t)  m  ,  (t)  -  {b  ,  (t)  +bx  (t)  }  m(  ^t) 

.  (B-6) 
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Appendix  C  (derivation  ) 


From  (3,2)  substituted  by  (3,3) 


V,  It) »  /o,  8;  it) 


kw/Zq,  yt) 

A 


<^t 


b/ci.\KW)+  iV«Xi*>  t  /q,uh |  ^ 

-/0„<Ait)  -J'Jk.ltij  Y,U)/;— ^  4WVtiiU) 


-JUk,1-*)  { X  W//q4  1-  i  WoHi  U ) 

_ <iW  v/q.?;W 

X,«)iY,lt)  ■>  /o.l^W  -t  tw) 


X.w//Q.t|,tt) 

-fc, 


S 


.  (C-1) 

Since  the  term 

_ I _  _  I 

XS*)  *  Y,  <*)  */50  {  Vt)  -t  V»  j  ^  +  4,^ 

X,U)  t  Y,(t) 

Q.  j  $,«)  -v  i  tt> }* 

a. 

|x,i-fc)  ■»  Yu>l 

<3^  {  3,  ^  +  "f>»  (-t)  | 


(C-2) 
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Substitute  (C-2)  into  (C-1)  and  take  Q--t<ao  #  then,  by 
same  logic  ; 


\  -  4  w^u,,wU,u) 

at  l^itj-v  4  )  d  l  *  h'  J  ‘ 

a 

dY'tt)=i37^7 

-  +  j  Yt-t)  dt 

_  atct)3,ft)lx,w  •*  Y,  w }  At 

I  J.tfc)  -t 

+  Ah^'X.U)J - - — -  - 4-\  «*vyA,  (.*) 

J  a'  1  3,W  •>  lud  -ft.  1  '' 

-AT  itiftt)  d  U/„  It)  ’ 

J  >  H. 


4Wus^) 

^ere  1  I?, 


The  ether  equations  were  derived  by  the  same  logic. 


the 
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g^o  v?.Q  too  ?c  o 


*  i  -  \ 


X-3CRLE=-5*  00E+00  UNITS 
V-SCRlE- I » 00E-03  UNITS 

MODEL-1 


Cl 


j 


N 


•  i 

HCH 

n~ 


US 


TIME 


*  l  —  £ 


X-SCfiLE«5.00E+00  UNITS  INCH. 
'■'-SCfiLS-5.00S.-04  UNITS  INCH. 

MODEL-1 

RUM  1  Q-2  US  TIME 


K-SCPLF.=5.00E+00  UNITS  INCH. 
y-3CPLF-l.00F.-03  UNITS  INCH. 


MODEL-1 


RUN  1 


M-l  US  TIME 


st> 


X-SCRLF.«5.00F.*00  UNITS  INCH. 
V- SCRLE-5.0CE-04  UNITS  INCH. 


M3DEL-1 


RUN 


h-2  US  TIM 


c 


sn 


go  cn 


M->2  iS  3e.ro  <mU 'cifti 


S;nC a.  ^  Qt.  <*r<  v 

T>v.  S  <^ra.?VKC<\l  0 

,’S  just  rfrindo W  flO\S€ 


CRLl=-5»  00E>00  UNITS  INCH. 

CRLC-'-'-L  OOE-iO  UNITS  INCH. 

VODEL-1 

RUN  1  M-12  N5 


*  2  - 1 


X-SCPILE>2.OOE>0O  units  inch. 

V-SCPLE-2. 00E-04  UNITS  INCH, 

MODEL-2 

RUN  i  Q-l  US  TIME 


X— 5CRLE--2 -  OCE+OO  UNITS  INCH, 
'v*-SCRLE=*l.  00E-04  UNITS  INCH, 


10DEL-2 


RUM 


n  o  lie  TIN  - 
1  J  ~  /  Uj  !  1  i  I  l_ 


go 


rt-SCRLF>2- OOE+OO  UNITS  INCH. 
y-SCRIl-'-.I .  00E-04  UNITS  INCH. 


MODEL-2 


M-l  US  TIME 


cro  nzo  °.2C  r?T0  oto 


K-3CRLE---2.  OOF.+OO  UNITS  INCH. 
V-2CRLE=5 •  QOE-05  UNITS  INCH. 

MODEL-2 


RUN  1 


k  A 


n 


1  'S 


u 


I  t  ,V1  _ 


*  z- 1~ 


K-SC,RLE>2.00E>00  UNITS  INCH. 

Y-3CRLF -2-  00F.-05  UNITS  INCH. 

MODEL-2 

PUN  I  M-12  US 


M  C 

i  ;  • 
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K-SCPLF-S. 00F>00  UNITS  INCH. 
'f-3CflLEI-2.00F.-04  UNITS  INCH. 


HODEL-3 


RUN  1 


G-l  US  TIME 


mo  l’to  ot:.  fcc»  ooo  roo  ?x 


N-SCRLE---5-  00E+00  UNITS  INCH, 
l'-3CRLr>2. OOE-04  UNITS  INCH, 


MODEL  ° 


■U 


RUN 


H-i  US  TIME 
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IUQ  2ir*  UdZ  a  lit,  QjO  (JjS 


— 

a 

1 

t 

i 

i 

1 

1 

1 

! 

j 

\ 

\ 

\ 

; 

j 

Old 


;:a 


020 


j25 


030 


«-3CfiLf«5.00E*00  UNITS 
y-SC, RLE-5.  00E-04  UNITS 

N0QEL-3 

RUN  1 


iCr  . 


MN-1 


I  |  o- 

Ub 


Th 


r 


X-5CfiL£"2 * OOE+OO  UNITS  INCH; 
'"'-3CRLE-2 •  GQE-04  UNITS  INCH. 

MODEL-3 

hv_  i 


Dl  !N  1 


IJS  TTMr£ 
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cn  cn 


CRLF.':2>  00EJ-00  UNITS  INCH. 

CRL: ~2*  ORC -04  UNITS  INCH. 

MODEL-3 

i  US  I  I  ME 


6* 


^  4-  ft  ~  L 


X-SCALE=2-  0OE-OO  UNITS  INCH. 
Y-SCALE-2- 00E-Q4  UNITS  INCH. 

M0DEL-4R 
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RUN  1 


C-i  US  TIME 


MO 


r-4 

a 


o 

j 


; 

i 


LT* 

o 


t\) 
L  i 


(  i 
O 

o 


’ 


K-SCRLF.--2 . OOc+OO  UNITS  INCH. 
V-SCRLEb2 • 00E-04  UNITS  INCH. 


MQQEL-4R 


RUN  1 


!_ 

t  I 


US  TIM: 
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OiO  0ir3  0*jQ  QdS 


iS-SCRLF.=-2. 00E*00  UNITS  INCH, 
V-SCRLE-5. 00F-04  UNITS  INCH, 


MGDEL-4R 


RUN  1 


MN-1  US  TINE 


7  l 


K-SCRLE>2.00E+00  UNITS 
''-SCflLE=l. 00E-03  UNITS 


MODEL-4! 


-1  o 


INCH. 

INCH. 


MV-1  OS  TIME 


-71 


c*tj  s?c  pxo  vrc:  r.To 


(  ) 


L  > 
_7 

t 


K-SCfiLD.2.00E*00  UNITS  INCH. 
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